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We discuss the bifurcation structure of homoclinic orbits in bimodal one dimensional maps. The 
universal structure of these bifurcations with singular bifurcation points and the web of bifurcation 
lines through the parameter space are described. The bifurcations depend on two parameters (codi- 
mension 2 bifurcations). We find the bifurcation lines exactly in a symbolic dynamics parameter 
plane and numerically in the parameter planes of a polynomial map and a piecewise linear map. 
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The aim of this article is to describe the structure of the bifurcations of some non-periodic orbits in a bimodal map. 
A more complete discussion of the bifurcations of orbits in bimodal, trimodal and four-modal maps will be given in 
ref. [1] (see also Ref. [2]). Bifurcations of periodic orbits in bimodal maps are investigated by several authors [3-7]. 

The bimodal map is a one dimensional continuous map Xt+i = /(^t ; a, h) where the function f[x) have one maximum 
point /(kci), and one minimum point f[xc^), and two parameters (a, 6) such that the two extremum points can be 

changed independently. We call x^^ and x^^ the critical points of the map. If x^^ < x^^ we denote the map -| h 

and if > x^^ we denote it — | — . In this letter we choose only to study the map -| h while the — | — map will 

be described elsewhere [1]. 

Periodic orbits are solutions of the equation 

ON 

2 3(a;) = /(™)(a;) - K = 0. 

^The solution of g{x) = is the topic of singularity theory [8] and the bifurcation of periodic orbits in the bimodal 
Omaps have codimension 2 structure, i.e. a universal form in a two-parameter plane. Non-periodic orbits are solutions 

|/"pf a slightly different equation and yields a different bifurcation picture but also with codimension 2. 

^ We choose the polynomial map 

Kt+i = x\ - axt^h (1) 

O 

^~^as an example of a smooth bimodal map, and as a continuous piecewise linear map we choose 

o 

O r 2a;t-Fa-F6-F2 if < -1 

. Kt+i = Ih- axt if -1 < < 1 (2) 

[ 2a;( - a -I- 6 - 2 if > 1 

^^These two very different bimodal maps illustrate which structures that are common for all bimodal maps, and the 
I'OViancfps of the structure when the smoothness oi f(x) is different. 
I We define symbolic dynamics [2,7,9-12] for the bimodal map St G {0, 1, 2} 

o 

^ r if < 

S( = < 1 if < < Kcj (3) 

[2 if > 

An orbit is identified by an infinite symbol string S = . . . S-2S-\SqS\S2 ■ ■ ■ and the points . . . , a;_2, JC-i, xq, xi, X2, ■ ■ ■ 
in an unstable orbit are uniquely defined by the symbol string. A point Xt of an orbit S has one symbolic future 
S(_|_iS(_|_2St+3 • • • and one symbolic past . . . St-2St-iSt. Since the bimodal map is not uniquely invertible a point Xt 
can belong to different orbits and have different symbolic pasts, but the point always has one unique symbolic future. 
We use the convention S1S2 • • • 5^ — {^1^2 • • • 

)°° to denote an infinite repetition of a finite length symbol string. We 
also use the notation . . . s„_i{s„, s^}s„_|_i . . . for the two symbol strings . . . s„_is„s„_|_i . . . and . . . s„_is^s„_|_i . . .. 

It is useful to define the following variable for the -|- — -|- bimodal map, following ideas of Milnor and 
Thurston [2,11,13]. 
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Pi 



= Si 

_ ( 1 if Si = or 2 
~ \ -1 if Sl = 1 



ist if pt-i = 1 

(2-.0 iip,.i = -1 (4) 
Pt-i if St = or 2 
-pt-i if St= 1 

CO 

t{xo) = 0.W1W2W3 • • • = ^ 

t=l 

where the base 3 real number < t < 1 is called the symbolic value. The kneading values are defined as 

Ki = t{Xc^) .gN 

Let Xt, t ^ ■[. . . , —2, —1, 0, 1,2,.. .} be the points of an orbit specified by the symbol string S and let the maximum 
and minimum values of the symbolic value be given by 

T™^^{S) = maxT(a;() with Xt < Xc^ 
T2™"(5) = mmT(a;() with Xt > x^, 

The orbit S is admissible if and only if 



^max 
'1 

K2 < t| 



,-min V / 



The values ki and K2 depend on the parameters a and h and we denote them the symbolic parameters of the 
bimodal map and the plane {ni,H2) the symbolic parameter plane. An orbit S is admissible in a rectangle in the 
(ki, K2) parameter plane. At the symbolic parameter values ki = 1.0 and K2 = 0.0 all orbits in the map exists. With 
these tools we can analyze the bifurcations of all orbits in the bimodal map. 

For simplicity we will here discuss only one family of homoclinic orbits. Other non-periodic orbits will have similar 
structures and these orbits can be analyzed in the same way. 

A homoclinic orbit in a bimodal map has a symbolic description 

Sa,b = 'ABA 

where A and B are finite symbol strings. We choose the simplest case with A = 1 

Si,B = IBl (8) 

with B = bib2 . . .bn- 

The bifurcation lines of Si^b with a few different strings B are drawn in the symbolic parameter plane (ki, K2) in 
Fig. 1. The two homoclinic orbits which exist for the largest area in the parameter plane are 121 and 12121 (in the 
notation introduced above T{1, 2}12T)_which both have rf = t(2T) = 0.2T = 5/6 and Tj™" = t(12T) = O.IOT = 7/18 
yielding the bifurcation lines K2 = 0.21 and K2 = 0.101. 

The homoclinic orbit is drawn in Fig. 2 (a) with a = 2.00794, b = 0.274 that are parameters close to the singular 
point in the parameter plane where both the critical points Xc^ and Xc^ belongs to the homoclinic orbit. In contrast 
to the bifurcation of periodic orbits is it important to notice which critical point that is first and which that is last 
along the orbit at the singular parameter point. In Fig. 2 (a) the orbit is first close to Xc^ then close to Xc^. There 
are four orbits bifurcating at this singular point with the symbolic description 

T{1,2}{0,1}2T (9) 

The structure of the bifurcation lines from the singularity depends on whether the line is associated with the first or 
the last critical point. The last critical point that the orbit visits yields one bifurcation line for all four orbits. In this 
example the line is ki = 0.21 in the symbolic parameter space. The critical point first visited by the orbit yields two 
difl'erent bifurcation lines; K2 = O.IOT for the orbits T{1, 2}12T and K2 = 0.02T = 5/18 for the two orbits T{1, 2}02T. 
These three lines are drawn in the symbolic parameter plane in Fig. 1 (a). The two K2 lines hit the ki line at two 
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different points, but these two points can be identified with one singular point because the line segment ki = 0.21, 
0.021 < K2 < 0.101 is not accessible. 

The bifurcation lines in Fig. 1 (a) are drawn in the (a, 6) plane in Fig. 3 (a) for the polynomial map (1). The 
main difference from the symbolic parameter plane is that the two bifurcation lines corresponding to the symbolic 
parameter lines K2 = 0.021 and K2 = 0.101 hit the line corresponding to ki = 0.21 approximately like the two branches 
of a parabola hit tangentially the end-points of a straight line. The area in the parameter space which gives a specific 
solution close to a singularity is often called a cusp because of the shape. In this case we have two different cusps; one 
which is a narrow sharp cusp and one which does not look like a cusp at all. The boundary lines of the second cusp 
have a common first derivative at the singularity but higher derivatives differ; the curve is but not C2. We believe 
that this is generic for smooth bimodal maps with /"(kcJ ^ 0. A simple way of understanding the cusp structures is 
to look at the solution of the simple equation 

Mfa{x)) = (10) 

where fa{x) and fb{x) are functions with one extremum value each. Choosing parabolas fa{x) = — a and fb{x) = 
x^ — b gives 

{x^ - a)2 - 6 = 

which has four real solutions for a > and < 6 < (a sharp cusp) and which yields at least two real solutions in 
the area bounded by 6 > when a > and b > when a < (a smooth cusp). 

The bifurcation lines for the piecewise linear map are drawn in Fig. 4 (a). Here the three bifurcation lines simply 
hit the singular point with different angles giving two cusps where the derivatives of the border lines are different 
at the singularity. The bifurcation structure of the piecewise linear map is related to the solution of (10) when the 
functions fa and are piecewise linear. Choosing fa{x) = |a;| — a and fb{x) = |a;| — 6 gives 

||a;| -a\-b = 

and this equation has four real solutions for < 6 < a and it has at least two real solutions in the area bounded by 
6 > when a > and b > —a when a < 0. These are two cusps of the same type as the cusps in Fig. 4 (a). 

Equation (10) captures the essence of a bifurcation where an orbit is born when this orbit does not have any 
point close to the critical points in an arbitrarily long future after visiting the two critical points once. The function 
fb{fa{x)) yields the same structure as obtained locally around an orbit visiting in any order: two times close to 
independent extremum points, and a finite number of times at points where the function is monotonically increasing 
or decreasing. An orbit is not close to any critical point (e.g. an unstable periodic orbit) can be considered to have 
constant points close to the singularity, and we choose a constant point in this orbit to be zero. The bifurcation of 
periodic orbits is different and is related to the solutions of an equation like 

fbifaix)) = X 

which gives different structures. The borderlines of a typical cusp for periodic orbits separates proportional to the 
distance from the cusp to the power 3/2 for smooth maps [8] compared with the narrower cusp of equation (10) with 
smooth functions where the distance between the borderlines increases as the distance from the singularity to the 
power 2. 

So far we have found that the structure of the singular point is different for periodic orbits and for homoclinic orbits. 
A further difference is that two singular bifurcation points for periodic orbits are always separated by a finite distance 
both in the symbolic parameter plane and in the parameter plane of a smooth map, while singular bifurcation points 
of homoclinic orbits may be arbitrarily close and may have bifurcation lines in common. 

In Fig. 1 (b) the bifurcation lines for the four homoclinic orbits 

T{0, 1}2{1,2}02T (11) 

are added in the symbolic parameter plane. These four orbits have the bifurcation line K2 = 0.021 in common with 
the orbits T{l,2}02Tbut have difl'erent ki bifurcation lines; kj = 0.2120T= 139/162 for the two orbits T{0, 1}2102T 
and Ki = 0.22021 = 149/162 for the two orbits 1{0, 1}22021. This yields a singular bifurcation point with three 
bifurcation lines but with ki and K2 switched compared to the structure of the bifurcation discussed above, because 
at this singular point the orbit (11) visits the point Xc^ before it visits Xc^. The homoclinic orbit close to the singular 
parameter value of the polynomial map is drawn in Fig. 2 (b) with a = 2.34865 and b = 0.26552. 

The structure of the bifurcation lines are drawn in the (a, 6) plane of the polynomial map (1) in Fig. 3 (b) and 
the structure of the singularity is the same as for the orbits (9). The bifurcation lines for the piecewise linear map 
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are drawn in Fig. 4 (b). This latter singular point is a finite distance away from the singularity discussed above but 
choosing for example the homoclinic orbits 1{0, l}2l(^™){l, 2}021 these will yield singular bifurcation points converging 
to the singular bifurcation point of the_orbit T{1, 2}{0, 1}2T since these have ki = 0.2l(2™+i)20T or ki = 0.2l(2™)202T 
which in the limit n ^ oo is ki = 0.21. These singular points will be arbitrarily close also in the (a, b) plane of the 
polynomial map and the piecewise linear map. The set of all singular bifurcation points on the line K2 = 0.021 is a 
complicated set of the Cantor type. Since each new bifurcation line in the Cantor set of lines again has new singular 
points in a Cantor set etc. this yields a complicated web of crossing bifurcation lines through the parameter plane 
connecting the homoclinic orbits of the form IBl. 

If the homo-/hetroclinic orbit ends in a periodic orbit which is created in a tangent bifurcation, for example a period 
3 orbit, then the cusp structure described above take place at the tangent bifurcation line in a tail of the swallowtail 
of the stable orbit. 

The parameter values when the critical point maps into an unstable orbit in the unimodal map is called Misiurewicz 
point and corresponds to the bifurcation lines described here. Some of these bifurcations corresponds to major changes 
in the structure of the attractor, such as band merging bifurcations and crisis bifurcations. To understand how the 
chaotic attractor changes with the parameters these bifurcation lines has to be described. 

In summary we have studied one example of the bifurcation structure for homoclinic orbits in bimodal maps. We 
find singular bifurcation points with cusps connected in a web of bifurcation lines through the parameter space. This 
is generic for all bimodal maps while the exact shape of the cusp depends on the smoothness of the function at 
the critical points. The bifurcation lines is found in a simple and exact way in the symbolic parameter space. The 
bifurcation structure for the homoclinic bifurcations in the two-dimensional Henon map is similar but not identical 
to what is obtained here and this is discussed in [14] 

The author is grateful for financial support from the Norwegian Research Council (NFR), and for interesting 
discussions with Predrag Cvitanovic, Jan Fr0yland and Tore M. Jonassen. 
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FIG. 1. Bifurcation lines for homoclinic orbits in the symbolic parameter plane. 



FIG. 2. Homoclinic orbits in the polynomial map (1) close to a singular bifurcation point, (a) 1{1, 2}{0, 1}21 for a = 2.00794 
b = 0.274, (b) I{0, 1}2{1, 2}02l for a = 2.34865 b = 0.26552. 



FIG. 3. Bifurcation lines in the parameter plane for the polynomial map (1). 



FIG. 4. Bifurcation lines in the parameter plane for the piecewise linear map (2). 
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